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Abstract 

We propose a new approach to the problem of the first passage time. Our method is applicable 
not only to the Wiener process but also to the non-Gaussian Levy flights or to more complicated 
stochastic processes whose distributions are stable. To show the usefulness of the method, we 
particularly focus on the first passage time problems in the truncated Levy flights (the so-called 
KoBoL processes), in which the arbitrarily large tail of the Levy distribution is cut off. We find 
that the asymptotic scaling law of the first passage time t distribution changes from i~( Q+1 )/ a -law 
(non-Gaussian Levy regime) to t~ 3 / 2 -law (Gaussian regime) at the crossover point. This result 
means that an ultra-slow convergence from the non-Gaussian Levy regime to the Gaussian regime 
is observed not only in the distribution of the real time step for the truncated Levy flight but also 
in the first passage time distribution of the flight. The nature of the crossover in the scaling laws 
and the scaling relation on the crossover point with respect to the effective cut-off length of the 
Levy distribution are discussed. 
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I. INTRODUCTION 



The first passage process or the first passage time (FPT) problem deals with the event 
where a diffusing particle or a random-walker firstly reaches a specific site at a specific time 
jl]. These FPT problems have been studied in various research fields, such as statistical 



hysics, chemistry [2] and biological neuroscience [31,14]. In finance, several authors 



flfl 



12] have analysed tick-by-tick data of the US dollar/ Japanese yen (USD/JPY) 



exchange rate and studied the FPT distribution for which the FPT is defined by the time 
that the rate firstly moves out from a given range. 



ii 



are reproduced 



12|]. That is, if 



Among these studies, the USD/JPY exchange rates of the Sony Bank 
from the market rates by using some rate windows with a width of 0.1 yen 
the USD/JPY market rate changes by more than 0.1 yen, the Sony Bank rate for USD/JPY 
is updated to the market rate otherwise it remains constant. In this sense, it is possible for us 
to say that the procedure for determining the USD/JPY exchange rate of the Sony Bank is 
essentially the first passage process. Despite many demands from various research fields and 
business in financial markets, one could obtain explicit analytical expressions or solutions of 
the FPT distribution only in very few cases. In addition, except for a few cases [2| , most of 
the analytical expressions are of the ordinary Wiener process (ordinary Brownian motion). 

Based on this fact, here we propose a new approach to the problem of the FPT or first 
passage processes. Our method is applicable not only to the Wiener process but also to 
the anomalous diffusion of the non-Gaussian Levy flights or more complicated stochastic 
processes. To show the usefulness of our approach, we particularly focus on the FPT prob- 



lems in the truncated Levy flights 
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221 ]. in which the arbitrarily large tail of the Levy 



distribution is cut off. Using the method, we find that the asymptotic scaling law of the 
FPT t distribution changes from a t~( a+1 ^ a -\&w (non-Gaussian Levy regime) to a t _3//2 -law 
(Gaussian regime) at some crossover point. This fact means that the crossover between 
non-Gaussian Levy and Gaussian regimes is observed not only in the distribution of the 
real-time step of the truncated Levy flight, which was reported by Mantegna and Stanley 



but also in the FPT distribution of the flight. Moreover, we give a scaling relation on 
the crossover point with respect to the effective cut-off length of the Levy distribution. The 
scaling relation enables us to predict the crossover point of the FPT distribution for a given 
truncated Levy flight. 
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This paper is organised as follows. In the next section, we explain general formalism of 
our method and apply it to the FPT problem for the Wiener process, for which the solution 
of the FPT distribution is well-known, in order to check the validity of our method. In 
Sec. Ill, we show that our method is widely useful for a class of stable stochastic processes. 
We derive the FPT distribution for Levy flight which includes Gaussian and Lorentzian 
stochastic processes as its special cases. For each stable stochastic process, we discuss the 
scaling law of the FPT distribution in the asymptotic regime. The analytical results are 
confirmed by computer simulations. In Sec. IV, we apply our method to the FPT problem of 
the truncated Levy flight and discuss the crossover in the scaling laws of the FPT distribution 
between non-Gaussian Levy and Gaussian regimes. The last section is a summary. 



II. GENERAL FORMALISM 



The problem we deal with in this paper is defined as follows. Let us consider the stochastic 
process : X k : k — 0, 1, • • • , T. For this time series, the FPT t is defined by t = minjfc > 
ko ; Xk = ±e}- Then, our problem is to obtain the distribution of t, namely, the first passage 
time distribution P(t). In other words, we evaluate the distribution of t, that is P(t), which 
is defined as the survival probability that the time series X k) starting from k = k keeps 
staying within the range [X ko — e, X ko + e] up to the time step k + t. The problem we 



are dea 
rate 
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ing with is motivated by the real mechanism of the Sony Bank foreign exchange 



121 ] . The Sony Bank rate is the foreign exchange rate that the Son y B ank offers 



with reference to the market rate. Basically, trades can be made on the web [13j while the 
market is open. The Sony Bank rate depends on the market rate but is independent of 
the customers' orders. If the USD/JPY market rate changes by e = ±0.1 yen or more, the 
Sony Bank rate for USD/JPY is updated to the market rate. For instance, for the stochastic 
process of the real market (what we call tick-by-tick data): Xq, X±, X2, ■ ■ ■ , Xt with < e 
and |X 2 | > e, the Sony Bank rate stays flat from the time k — to k — 1, as the market 
rate is in the range of e = ±0.1 yen based on the market rate at k = 0. When the market 
rate exceeds the range of e = 0.1 yen at k — 2, the Sony Bank rate is updated to the market 
rate. Obviously, the time interval t = 2 — = 2 here corresponds to the FPT we explained 
above and it is worth while for us to evaluate its distribution P(t) to investigate statistical 
property of the Sony Bank USD / JPY rates. 
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To calculate the FPT distribution P(t) for the time series X k , we define the probability 
P(l), which means the probability of the FPT is t — 1 as 

T 



1 T 

P(l) = lim - ©(l^+i - M - e) (1) 



fc=i 



where G(x) means the Heviside step function, namely, G(x) = 1 for x > and 0(x) = for 



x < 0. We usually solve a 



dnd of 



ractal) Fokker-Plank equations under some appropriate 



boundary conditions 0, 3], ^, 14,[l5,Q or use the so-called image method Q, Q to discuss 
the FPT problem. However, as we saw in equation ([I]), our approach is completely different 
from such standard treatments. To evaluate the FPT (probability) distribution, say -P(l), 
we directly count the number of t = 1, namely, Mi = J2k=i ©(l^fc+i — Xt| — e) appearing 
within quite a long time interval T. We might choose T as a time interval during which 
the market is open. Then, the ratio Ni/T should be expected to converge to -P(l) as T 
tends to infinity. This is the meaning of equation and is our basic idea for evaluating 
the FPT distribution. From our method, to evaluate the FPT distribution by counting Mt 
(t = 1, 2, ■ ■ ■ ), the probability -P(2) is also given by lim T ^ 00 (A2/T), that is to say, 

T 



P(2) = liui i J2e(\X k+2 -X k \-e)-P(l) 



k=l 

T „ T 



= hm Lj2e(\X k+2 -X k \-e)-lim 1 £ G(\X k+1 - X fc | - c). (2) 

k=l k=l 

In the same way as the probability -P(2), the probability P(3) is obtained as 

1 T 

P(3) = lim ^6(|X fe+3 -X fc |- e )-P(l)-P(2) 



fe=l 

T T 



= lim ±J2e(\X k+3 -X k \-e)- lim ± £ 0(|X fc+2 - X fc | - e). (3) 

k=l k=l 

We should notice that the probability P(l) was cancelled in this expression Q. Thus, we 
easily generalise this kind of calculations to evaluate the distribution P(t) by repeating the 
above procedure as follows. 

T T 

P(t) = lim ± d Xfc +' " Xk \ " e ) " T lim t E d X ^i " M ~ e ) ( 4 ) 

fe=l fc=l 

where P(l), • • • , P(i — 2) were all cancelled in this final formula pi) . This equation (j3J) is 
the starting point of our evaluation. At a glance, this equation seems to be just a definition 



of the FPT distribution; however, for some classes of stochastic processes, we can derive the 
explicit form of the FPT distribution from this simple equation. In the next subsection, we 
derive the FPT distribution for the Wiener process as a simple test of our method. We stress 
that our approach helps as an intuitive account for the first passage process and derivation 
of its distribution. 



A. A simple test of the method for Wiener stochastic processes 

To show the validity and usefulness of our method, we derive the FPT distribution from 
the above expression (j3J) for Wiener stochastic processes (Brownian motion). The ordinary 
Wiener process is described by Xt+i = X t + Y t , where the noise term Y t obeys the white 
Gaussian with zero- mean and variance a 2 . Then, we should notice that the difference 
St = X k+t — X k is rewritten in terms of the sum of the noise terms Y t as S t = X^=o^fc+j- 
As is well-known, as the Gaussian process is stable, St obeys the Gaussian with zero-mean 
and ((S t ) 2 ) = to 2 variance. Using the same argument as S t , St-i also obeys the Gaussian 
with zero-mean and ((S t -i) 2 ) = (t — l)a 2 variance. Therefore, the FPT distribution P(t) 
derived by equation leads to 

P(t) = Q(p:S t ,S t ^) (5) 

POO POO 

n(p:l,m)= / p(l)Q(\l\ - e)dl - p(m)Q(\m\ - e)dm (6) 



when we assume that the underlying stochastic process is ergodic, namely, that time av- 
erage and ensemble average coincide. For the ordinary Wiener process, the probability 
distributions for St and St-i are Gaussians with zero mean and variances ta 2 , (t — 1)<j 2 , 
respectively. Thus, we easily evaluate the integral appearing in (jSJ) after substituting 
p = PG (S t ) = (l/V2^ 2 i)e- s ?/ 2 ° 2t ,p G (S t - 1 ) = (l/ v /2vra 2 (t-l))e- 5 -i/ 2CT2 C- 1 ) and obtain 
P(t) = n(p G :S t ,S t - 1 ) as 

™- a {*fe)-*(vr*)} 

where we defined a = e/er. The function H(x) is defined by H(x) = dze~ z2/2 /V2^. We 
should keep in mind that the above result is valid for discrete time t, however; it is easy for 
us to obtain its continuous time version by replacing t — > t, t — 1 — > i — At and evaluating 



5 



(J7J) in the limit of At — > 0. Then, we have 

P it) M = 2ff (A) _ 2ff (_^_) = 2 |* (_l) Ai + 0((At) ,, (8) 

Thus, the FPT distribution for the ordinary Wiener process in the continuous time limit is 
given by 

P(t) = 2i-H ( 4= J = ■ (9) 



This well-known form is expected the inverse Gaussian distribution [18( for the FPT dis- 
tribution of the ordinary Wiener process and is often observed in the so-called inter-spike 
interval (ISI) of the integrate-and-fire model for neural networks ay 

. Therefore, in the 

asymptotic regime t — > oo, the FPT distribution for the Wiener process obeys t~ 3 / 2 -scaling 
law. From the above discussion, we found that our new approach based on direct counting 
of the FPT to obtain the FPT distribution is effective and gives a well-known solution for 
the ordinary Wiener process. 

Before we move to the main section, we should comment on the much shorter derivation 
of the above formulation ([9]). Let us define our stochastic process by x(t) and the probability 
density function finding the process in x at time t provided that it was in x at time by 
p(x ,t\x, 0). Here, we assume that p(x ,t\x, 0) is symmetric with respect to x = 0. We also 
define the FPT with absorbing barriers x = ±e by r. Then, the complementary cumulative 
distribution of r is given by 

W{x,t) = P r (r>t) = J p(x ,t\x,0)dx (10) 

and the cumulative distribution function leads to 

/oo 
p(x ,t\x,0)dx . (11) 

Thus, the FPT distribution is obtained by 

P(x ,t) = (12, 



This equation coincides with equatio n ([9 



derivations are found in the references 



24 



We should also mentioned that the similar 



25|. 
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III. STABLE PROCESSES AND THEIR FPT DISTRIBUTIONS 



We stress that our method is widely applicable to stochastic processes whose distributions 
are stable. Stable processes are specified as follows. If stochastic variables Yi {% — 1, • • • , N) 
are identically independently distributed from p(Yj), the Fourier transform of the sum of the 
Yi, namely, S n = Y,jLi Y j is g iven by 

Ml) = WW (!3) 

where <f{q) is the Fourier transform of the stochastic variable Y, namely, the characteristic 
function and defined by 

/oo 
p(Y) e iqY dY. (14) 
-oo 

Then, the stochastic process Yi is referred to as a stable process. Strictly speaking, equation 
(I13p is a possible definition of infinitely divisible random variables and not of stable random 
variables. A stable random variable is infinitely divisible and stability refers to the invariance 
of the distribution with respect to convolutions. 

It is obvious that for this class of stable processes, the FPT distribution is easily obtained 
by our method because the probability distributions p(S t ) and p(S t -i) to evaluate Q(p : 
S t , S t -i) in ([5]) are defined explicitly. In the next subsections, we show several results from 
our new approach. 



A. Lorentzian stochastic processes 

As a first simple example of the stable distributions, let us think about Lorentzian stochas- 
tic processes: X i+1 = X t + Y t , where the noise term Y t obeys the following white Lorentzian 



p(X*) = i-r-v-2- ( 15 ) 



7 1 

7T 7 2 + Y t 

Then, the characteristic function of the stochastic variable Y t is given by 



poo iqYt 

^^lL^TV dYt = e (16) 

By using the convolution of the Fourier transform for the variable S n = ^12=1 we have 
tPnil) — [y 9 (Q , )] ri — e _n7 ' 9 '. Therefore, the inverse Fourier transform of (p n (q) leads to the 



probability distribution of the sum of noise term as follows. 

By substituting this probability PLorentz(S t ) , PLorentz(S t -i) into equations (j5J) and we 
obtain the FPT distribution for the Lorentzian stochastic processes as follows: 



00 



P(t) — Q(pLorentz '■ S t , S t -l) — 2 / — „ — 2 ^ 



* (7*) 2 + s? y« t ( 7 (( - + s|_, 

2 tan- 1 f -i-) --tair 1 (-) (18) 



7T \t— 1/ 7T \t y 

where we defined 6 = 6/7. This is a result for discrete time steps; however, its continuous 
version is easily obtained by using d{ian^ 1 {y))/ dy = (1 + y 2 )" 1 as follows. 

9 (2 .mi 2b 1 

From this result (TT9T) . we find that the FPT distribution for the Lorentzian stochastic pro- 
cesses obeys Lorentzian. In the asymptotic regime t — > 00, the FPT distribution for the 
Lorentzian stochastic processes obeys the t _2 -scaling law. 

B. Anomalous diffusion of Levy flight 

We next consider the case of Levy stochastic processes whose noise term Y t of the stochas- 
tic process = X t + Y t obeys the following Levy distribution : 

1 r°° 

PL6vy(Y t ) = - / e-^' Q cos(qY t ) dq (20) 



We should keep in mind that the above distribution (1201) is reduced to the Wiener stochastic 
process (a = 2) and the Lorentzian stochastic process (a = 1) as its special cases. As this 
process Y t (t — 1, • • • , n) is also stable, the sum of the noise term S n = YTj=i Yj nas the 
following probability distribution 

PLevy(S n ) = - / cos(qS n ) dq. (21) 

TT Jo 

Now, we can derive the FPT distribution by substituting PLevy(St) , PL6vy(St-i) into equations 
(ED and (ED as 



P(t) = &{pLevy, S t , S t -l) 

cy POO POO Q /'OO /*00 

= - / / e^l" cos(qS) dq-- dS e^-VM* cos (qS) dq. (22) 

* Je J0 TT J e Jo 



Expression of a continuous time version ([22"]) is obtained from the derivative of the above 
discrete time distribution P(t) with respect to t as 

/"OO /*00 *\ r\ /*oo /•OO 

dS e-^l Q cos(qS) dq } = — ^ dS \q\ a cos(qS) dq. 

Je J0 J 71 Je JO 

(23) 

In the asymptotic regime t — ► oo, by replacing the variable as tq a = Q and after some simple 
algebra, we obtain 



P(t) = ±25 / dS / Q 1 /a e -^a,s 



7TO; 

-(a+l)/a 



= ^(a, 7 )r (a+1)/Q (24) 



7ra 2 ^ ' (2/)! V « 

where the r(x) means Gamma function and p is a constant of order 1. Then, we should 
notice that the above scaling law is consistent with both the Wiener stochastic process 
(P(t) ~ t~ 3 / 2 for a = 2) and with the Lorentzian stochastic process (P(t) ~ t -2 for a = 1) 
we discussed in the previous subsections. 

As we saw in the above stochastic processes, our new approach based on direct counting of 
the FPT is widely useful for the class of stable stochastic processes and the final expressions 
containing at most just two integrals. Moreover, our approach can be applied to the 
FPT problems with a surprisingly wide variety of absorbing boundary conditions. This is 
one of the advantages of our method over other approaches based on analysis of Fokker- 
Plank equations. To show the advantage, in the next section, we apply our method to much 
more complicated stochastic stable process. 



C. Computer simulations 



To check the theoretical prediction of the power-law exponent for the FPT distribution, 
we perform computer simulations. Then, to generate the additive noise Y t for each time 
step, the Levy stable distribution p(Y) is needed. As shown by Umeno |26[], the Levy stable 
distribution is obtained by a superposition of a chaotic map such as 



n+l 



V \ a 

\ 1 n\ 



1/a 



sgn(y„- — ) (0<a<2), 

n 



(26) 
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namely, for independently selected initial conditions i = 1,---,N, the quantity 
YliLi /N l l a obeys a Levy stable law with a, that is, 

1 f°° 

p(Y) = - / e~ qa cos(qY)dq. (27) 

In FIG. [TJ we plot the Levy distribution with a = 1 (Lorentzian) and a = 1.2 obtained 
by the superposition of the chaotic map (T26]) . We set the number of the superposition of 
the chaotic map N = 10000. From this figure, we find that the distributions obtained by 




-10 -5 N/ 5 10 -10 -5 N/ 5 10 



FIG. 1: Levy distributions obtained by the chaotic map ([26| with a = 1 (Lorentzian: left) and a = 1.2 
(right). The broken lines are corresponding theoretical results given by (|27|) . We set the number of the 
superposition N = 10000. 

the simulations are in good agreement with the corresponding analytical expressions (j27p . 
Keeping these results in mind, we use the sampling point from the superposition of the 
chaotic map (T26]) as the additive noise Y t in the stochastic process for each time step. Then, 
we should notice that one should choose a large value for the width of the rate window e, 
say, e ~ 100 to investigate the tail of the FPT distribution. This is because it is a rare event 
for the random walker to stay in the range [— e, e] and it takes quite long computational 
time for us to obtain the tail exponent (relatively long FPT) if e is small. Here, we choose 
e = 120 and evaluate the power-law exponent of the FPT distribution for the Levy processes 
with a = 1 and a = 1.2. The results are shown in FIG. [2J From this figure, we find that 
the power-law tails have almost the same exponents as those predicted by our theory. Of 
course, the tail region is too noisy to conclude that the simulations are completely consistent 
with the theory. However, at least, we may say that they are not in disagreement. Thus, 
the computer simulations provided us with a justification of our theoretical formulation. 

10 




t 



a=1.2 




t 



FIG. 2: The tail exponents of the FPT distributions for a = 1 (left) and a — 1.2 (right). The values of the 
power-law exponent predicted by our theory (the slope of the broken line in each panel) are —2 for a = 1 
and —1.8333 for a = 1.2, respectively. We set the width of the rate window e = 120. 

IV. CROSSOVER IN SCALING LAWS OF FPT DISTRIBUTIONS 



In the previous section, we showed our new formulation is effective and much more simpler 
than the approach of the (fractal) Fokker-Plank equations [14| to obtain the FPT distribution 
for stable stochastic processes. We actually found that the FPT distribution of the general 
non-Gaussian Levy stochastic process specified by parameter a is obtained and its scaling 
behaviour in the asymptotic regime t — > oo is t~( a+1 )/ a -iaw. In this section, we show that 
our formalism is also useful in obtaining the FPT distribution for the so-called truncated 
Levy flight (the so-called KoBoL processes from Koponen, Boyarchenko and Levendorskii 



22 



27|, 



281]). for which it is well-known that the crossover between the a Levy and a Gaussian 
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1231 ]. In this section, 



regime in the distribution of the real time step takes place QJ 
we show, by using our method based on direct counting of the FPT, this kind of crossover 
in scaling laws is also observed in the FPT t. 



The characteristic function for the truncated Levy flight is defined by 

(A 2 + |g| 2 )"/ 2 cos (a tan- 1 (M)) _ A a 



T(q) 



exp 



-7- 



cos 



\ 2 ) 



(28) 



where tan 1 (|g|/A) ~ 7r/2asA — ► 29(. Therefore, in this limit A — >• 0, the above equation 



( |28|) is reduced to T(q) = e 7 '' ? I Q . Obviously, this is identical to the characteristic function 
of the conventional Levy flight as we already saw in the previous section. For this reason, 
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a non-zero value of A controls the cut-off width of the truncated Levy flight. We should 
notice that one could also use a hard cut-off version of the truncation scheme [^(J, namely, 



PTLFiXt) = p^yiXt)^- 1 - \Y t \). 



(29) 



However, for its mathematical simplicity, we use the soft cut-off version of the truncation 
scheme which will be explained below. 

By the convolution of the Fourier transform, we can show 



T n (q) = [T(q)} n = exp 



-717- 



(A 2 + \q\ 2 ) a/2 cos (a tan- 1 ( M ) ) - A" 



cos 



\ 2 ) 



(30) 



and then the sum of the noise term S n = Y^Jj=\ Yj °f the truncated Levy flight obeys the 
following probability distribution: 



1 f°° 

PTLF(S n ) = — / dqe-^[T(q)Y 



TV 



dq cos(qS„ 



(A 2 + Ig^^cosfatan- 1 (M) ) - A c 



-wy- 



cos 



( 1Ta \ 
\ 2 / 



(31) 



Substituting these probability distributions PTLF(S t ) , pTLF(S t -i) into equations and ([5]), 
and taking the derivative of P(t) with respect to t, we obtain the FPT distribution of the 
truncated Levy flight for continuous time case as 



7T 



x exp 



dS I dq 
'0 



(A 2 + |g| 2 ) Q / 2 cos(atan- 1 (f 



-jt- 



(A 2 + |g| 2 ) Q/2 cos (a tan" 1 



cos 



I 7TCJ \ 

\ 2 ) 



A c 



cos 



\ 2 ) 



cos(qS). 



(32) 



Thus far, it has been difficult to perform the above two integrals with respect to S and q 
analytically to obtain a compact form of the FPT distribution. However, numerical integra- 
tions of equation (132]) enable us to proceed to it. In FIG. [3J we show the scaling plot of the 
FPT distribution for the truncated Levy flight with a = 1.2 for several values of A. From 
these three panels in FIG. IU we find that the scaling law of the FPT distribution changes 
from = £-(i-2+i)/i.2 _ j-1.83 to r 3/2 at gome crossover points t x ~ 44 (A = 0.5), 19 

(A = 1.0) and t x ~ 8 (A = 2.0). To obtain useful information about the crossover point t x , 
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a=1.2 A=1.0 
0.33 x 1 83333 
0.115 x- 1 5 





a=1.2 A=2.0 * 




0.30 x- 1 83333 




0.15 x- 1 5 
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FIG. 3: Scaling laws of the FPT distribution for the truncated Levy flight. We set a = 1.2 and A = 0.5 
(upper left panel), 1.0 (upper right panel) and A = 2.0 (lower panel). We find a clear crossover between Levy 
and Gaussian regimes. The scaling laws change at the crossover point t x ~ 44 (A = 0.5), 19 (A = 1.0) and 
t x ~ 8 (A = 2.0). Below the crossover point t x , the scaling laws are those of the Levy : f~( a + 1 )/ a — £~ 183 5 
whereas, above t x , the scaling laws become those of the Gaussian : t~ 3 / 2 . In these three panels, we find 
that the scaling relation : t x (A = 0.5)/i x (A = 1.0) = t x (A = 1.0)/f x (A = 2.0) = 2 a ~ 2.297 holds for the 
crossover points. 



we evaluate the asymptotic form of the FPT distribution (1321) for both \q\ < A (Gaussian 
regime) and |g| > A (non-Gaussian Levy regime). 

For Gaussian regime \q\ < A, by replacing the variable q with Q as 
7to 2 A Q - 2 g 2 /(2cos(vra/2)) = Q, that is, 



'2Qcos(f) tl/2 



7« 2 A 



2 Ao-2 



(33) 
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we obtain 



P(t) = V(a,7,A)r 3/2 

y/2 (ja 2 A a ^ 
i> a, 7, A = — 

7T V COS (^) 



-1/2 



g 1/2 e Q cos 



7a 2 A 



2 Aa-2 



2cos(^) 



Q 1/2 S 



(34) 
(35) 



It should be noted that this t 3,/2 -law is valid for \q\ < A. From the equation (1551) . this 
condition reads 



2Q cos(^) 



(36) 



7« 2 A a 

On the other hand, for \q\ > A, that is to say, for t < t x , the FPT distribution ( 1521) is 
evaluated as 

Vjcos(^) 



P(t) ~ 

7T 



7T 



diS 1 / 
'0 



cos 



( na \ 
\ 2 ! 



cos(qS) 



dS 



dq \q 



« -7kl Q * 



cos(qS). 



(37) 



P(t) 



This result is identical to the FPT distribution for the conventional Levy flight, which is 
defined by equations (T24l and ()25p . and was already obtained in the previous section. 

Let us summarise the result for the scaling laws of the FPT distribution for the truncated 
Levy flight. 

f-( a + 1 )/ a (f < 1 ■ non-Gaussian Levy regime) 

(38) 

t -3 / 2 (t >t x '■ Gaussian regime) 
We should bear in mind that the crossover point t x obtained by (1361) contains integral vari- 
able Q. Therefore, it is hard to say that t x is well-defined. To delete the Q-dependence of t x , 
we consider the ratio of t x (A) and t x (2A). From equation (1361) . we obtain £ x (A)/t x (2A) = 
2 a , namely, t x (A) = 2 a t x (2A). Let us check this scaling relation for the result we obtained 
in FIG. El For a = 1.2, the relation reads t x (A) = 2 L2 t x (2A) = 2.297t x (2A). This relation 
predicts the crossover point t x (A = 0.5)/t x (A = 1.0) = t x (A = 1.0)/t x (A = 2.0) = 2.297, 
which is very close to the results obtained in FIG. [51 namely, 44/19 ~ 2.316 and 19/8 ~ 
2.375. The small difference is probably because of impreciseness of numerical integrations 
appearing in equation (1521) . 

The relation t x (A) = 2 a t x (2A) for successive values of A and 2A is easily extended for 
the relation between A and 5 A (5 > 1) as follows. 

t x (A) = 5 a t x (5A) (39) 



14 



This scaling relation for the crossover point t x in the scaling laws of the FPT distribution of 
the truncated Levy flight t x is one of the main results in this paper. From this result (1391) . 
we find that the crossover point t x increases rapidly as the effective cut-off length I = (<5A) _1 
also increases as 

*x(0 = (r) t x (l ) (40) 



Jo, 

where we set l = A -1 . Therefore, we conclude that the crossover between non-Gaussian 
Levy and Gaussian regimes is observed not only in the distribution of the real-time n flight, 
which was reported by Mantegna and Stanley 20J, but also in the FPT distribution of the 
truncated Levy flight. 

In the study by Mantegna and Stanley 20], they investigated the stochastic variable 
z n = Y^k=i x ki where x k = X k — X k _ x obeys the truncated Levy flight. They evaluated the 
probability of return P[z n = 0) and found that the P{z n = 0) obeys the Gaussian n _1//2 - 
law in the large real time step n regime. In this section, it was shown that this ultra-slow 
convergence from Levy regime to the Gaussian regime is conserved even if we consider the 
first passage process of the truncated Levy flight. The relation between their results and 
ours is clearly understood as follows. 

For a given time interval t of the first passage process of the truncated Levy flight, the 
time series of the variable x k = X ko+k — X ko+k _i behaves as x\, x 2 , ■ ■ • , x t , where k is an 
origin for the measurement of the interval t. Then, from the observation by Mantegna and 
Stanley, the sum z t = $^t=i x k = X ko+t — X ko obeys a Gaussian with zero- mean and variance 
t if the time interval t is large enough, that is, t > t x . Then, the probability of return is 
given by P(zt = 0) ~ t^ 1 ^ 2 . In other words, for t > t x , it takes quite long time for a 
random- walker to escape from the region [X ko — e,X ko + e], and the time t for escaping 
guarantees that the central limit theorem works to make the variable z t a Gaussian. As a 
result, the FPT distribution P(t) should follow the corresponding Gaussian t~ 3 / 2 -law from 
our argument for the case of the Wiener process (Q. On the other hand, if the interval t is 
smaller than the crossover point t x , the central limit theorem for z t does not work and z t is 
no longer a Gaussian. Then, as we checked, the FPT distribution P(i) obeys t~^ a+1 ^ a -\a,w 
of the Levy flight. 
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V. SUMMARY 



In this paper, we proposed a new approach to evaluate the FPT distribution. Our method 
is based on direct counting of the FPT. We show that our approach gives an explicit form 
of the FPT distribution for stable stochastic processes. Actually, for Wiener (Brownian 
motion), Lorentzian and Levy stochastic processes, our method was demonstrated. Thanks 
to the mathematical simplicity of our method, it becomes easy to grasp the intuitive meaning 
of the FPT distribution and to tackle more complicated stochastic processes. As an example, 
we discussed the FPT distribution of the truncated Levy flight (the KoBoL process). We 
found a clear crossover between non-Gaussian Levy and Gaussian regimes in the scaling 
laws of the FPT distribution. We found the scaling relation on the crossover point t x with 
respect to the effective length I of the cut-off as t x (l) = (l/lo) a t x (lo) with l = A -1 . 

Very recently, Koren et.al. [30] investigated not only the FPT distribution but also the 
first passage leapover (FPL) distribution under a single absorbing boundary condition. A 
relatively new concept, the FPL is defined as the flight length for a random walker to move 
beyond the single boundary (a target). Our system in this paper possesses two boundaries 
(in this sense, our process might be referred to as a first exit process) ; however, it might be 
possible to apply our analysis to the problem in order to discuss the FPL distribution. This 
will be addressed in future work. 

We hope that beyond the present analysis for the Sony Bank rate, our approach might be 
widely used in many scientific research fields, especially in the field of econophysics including 
financial data analysis. 
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